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We examine the noise-induced transition in a fluctuating bistable potential of a driven quantum 
system in thermal equilibrium. Making use of a Wigner canonical thermal distribution for descrip- 
tion of the statistical properties of the thermal bath, we explore the generic effects of quantization 
like vacuum field fluctuation and tunneling in the characteristic stationary probability distribu- 
tion functions undergoing transition from unimodal to bimodal nature and in signal-to-noise ratio 
characterizing the co-operative effect among the noise processes and the weak periodic signal. 
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I. INTRODUCTION 

The role of the noise in nonlinear dissipative systems has been an intriguing issue over the last two 
decadesi*2iiii£i§iL£i2ii2iiiiI2iiiiI2ii£. Because of their potential applications, various noise-induced processes in zero 
dimensional and spatially extended systems, such as, stochastic resonance^Sii, coherence resonance^, thermal 
ratchelpiii&iii, noise-induced front propagation 8 , noise- induced patten format ion&iSiiAiiSii&ii etc have been the sub- 
ject of intense investigation in widely different areas of physical, chemical and biological sciences. An important early 
endeavour in this direction is the noise-induced transitions^ which plays a key role in understanding of the transition 
of unimodular character of stationary probability distribution to bimodal one as one varies the strength of external 
multiplicative noise in a bistable potentialS2S2ii. The theory has been further advanced^S^ in several directions to 
include the aspects of localization, the effect of correlation between the noise processes, the cooperativity between the 
noise processes and the signal and the associated transient characteristics. For example, it had been shown earliest 
that a system in a bistable fluctuating potential driven by a very slow modulation or static tilting in presence of an 
additive noise exhibits localization in one of the wells. The stationary probability distribution is markedly sensitive 
to the variation of the correlation between the additive and multiplicative noises. Keeping in view of these obser- 
vations on noise-induced transition in classical systems it is imperative that quantization2£L2ii22i2i is likely to affect 
these features since tunneling and other generic quantum effects start playing significant role in a wider context. As 
prototypical experimental systems^ it is now possible possible to confine electrons in semiconductors, for example, in 
a quantum dot coupled to an environment via point contacts using nanolithographic techniques to explore in detail 
the quantum transport in relation to localization and tunneling. 

The object of the present paper is to study noise-induced transition in a fluctuating bistable potential in a quantum 
system driven by a sinusoidal slowly varying field in presence of additive thermal noise of the heat bath. Our aim 
here is to explore the characteristics of the stationary probability distribution functions to understand the nature of 
noise-induced transition in presence of generic quantum effects like tunneling and vacuum fluctuations. In addition 
we also look for how these changes in the stationary probability distribution functions are reflected in the nonlinear 
co-operative effect between the noises and the signal. 

II. QUANTUM STOCHASTIC DYNAMICS OF A DRIVEN SYSTEM 

A. General aspects 

To derive quantum Langevin equation from a microscopic picture we consider the well-known standard system- 
reservoir model with following form of Hamiltonian^i 
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Here x and p are the coordinate and momentum operators of the particle and {qj,Pj} are the set of coordinate and 
momentum operators for the reservoir oscillators coupled linearly through the coupling constants Kj(j = 1,2,...). The 
potential V(x, t) is due to the external force field for the Brownian particle. The coordinate and momentum operators 
follow the usual commutation rules {x,p} = ih and {<ji,Pj} — ihSij. Eliminating the bath degrees of freedom in the 
usual way we obtain the operator Langcvin equation for the particle 

mx+ [ dt'-Y(t-t')x(t') + V'(x,t) = F(t) (2) 
Jo 

where the noise operator F(t) and the memory kernel j(t) are given by 

F(t) = \{Qj(0) - £(0)}«j cosojjt + «] /2 pj(0) sinu jt (3) 

3 

and 

l(t) =^/SjCOSWjt (4) 

3 

respectively, with kj = uq. Very recently it has been shown [28-30] that on the basis of quantum mechanical average 
(...) over the bath modes with coherent states and the system with an arbitrary state Eq.(2.2) can be cast into the 
form of the generalized quantum Langevin equation. 



mi + f dt'j(t - t')x(t') + V'(x, t) = f(t) + Q(x, (Sx n )) (5) 
Jo 

where the quantum mechanical mean value of the position operator (x) — x. Here the quantum dispersion term Q 
to the potential, is given by 

Q(x,(6x n )) = V'(x,t)-(V'(x,t)) (6) 
which by expressing x{t) = x{t) + Sx(t) in V(x, t) and using a Taylor series expansion around x may be rewritten 



as 



Q{x,{5x n )) = -Y J - ] V n+1 (x,t)(8x n ) (7) 



i>2 



The calculation of Q rests on the quantum correction terms (Sx n ) which one determines by solving a set of quantum 
correction equations (as discussed in the later part of this section). Furthermore the c- number Langevin force is given 
by 

/(*) = E [{&(°)> - <*(°)>}*j cos ^ + «} /2 fc(°) sin ^*] (8) 

3 

which must satisfy noise characteristics of the bath at equilibrium , 



(f(t))s = (9) 



(f(t)f(t')) s = ^ f coth ) cos ^(* - ( 10 ) 
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Ea. i|10|) expresses the quantum fluctuation-dissipation relation. The above conditions (2.9-2.10) can be fulfilled 
provided the initial shifted co-ordinates {(<Zj(0)) — (5(0))} and momenta (pj(0)) of the bath oscillators are distributed 
according to a canonical thermal Wigner distribution 3 ^*^ of the form 

pm,m - «o))], fe (o))) - *» P {- ife<°>>-+w>>-^>>n 

so that the statistical averages (...) s over the quantum mechanical mean value O of the bath variables are defined 

as 

(Oj) s = J O j P j dfo(0)) d{&(0)> - (x(0))} (12) 

Here n(u>) is given by Bose-Einstein distributions (e^# — Pj is the exact solution of Wigner equation for 

harmonic oscillatoi^S^ and forms the basis for description of the quantum noise characteristics of the bath kept in 
thermal equilibrium at temperature T. N is the normalization constant. In the continuum limit the fluctuation- 
dissipation relation l|10|) can be written as 

1 fill} 

(f(t)f(t'))s = 2 J q k(cj) p(cu) hw coth(— ) cosu;(t - t>) (13) 

where we have introduced the density of the modes p(u>). Since we are interested in the Markovian limit in the 
present context, we assume k(uj)p(uj) = —7, and Eq.O then yields2i 



(f(t)f(t')) s = 2D q S(t - f) (14) 



with 



1 fkjQ 

D q = -jhw coth — (15) 



u>o refers to static frequency limit. Furthermore from Eq.Q in the continuum limit we have 

1 (t-t')= 1 5(t-t>) (16) 
7 is the dissipation constant in the Markovian limit. In this limit Eq.JSJ) therefore reduces to 

mx + 7 i + V'(x, t) = f(t) + Q(x, (6x n )) (17) 

In order to consider the stochastic dynamics in a bistable system with fluctuating barrier height we consider the 
potential as, 
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, 2 + -/ (18) 



V(x,t) = 

£(t) is the external fluctuations with zero mean and is characterized by the following equations 



ffl = ; WW) = 2Q/«(* - 1') (19) 

where Qi is the strength of noise. Care must be taken to distinguish between the quantum mechanical mean (O) (= O), 
the statistical average over the quantum mechanical mean (0} s denoting thermal bath average and the statistical 
average over the external noise O. 

Attention is restricted here to overdamped Langevin equation in one variable x so that we write 



7± = (ax - 6a; 3 ) + x£(t) + f{t) + Q(x, (Sx n )) 



(20) 
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Q(x, (8x n )) being the correction term due to nonlinearity of the potential is too small to destroy bistablity or the 
symmetry of the system. If we now apply an weak periodic forcing to the particle, the double-well potential is tilted 
asymmetrically up and down and the corresponding stochastic differential equation reads as 



jx = (ax - bx 3 ) + x£(t) + f(t) + Q(x, (Sx n )) + Acosujt 
The Fokker-Planck equation corresponding to Eq. 1)21(1 is given by 



(21) 



dP(x,t) _ d_ 
dx 



where 



and, 



B(x,t) 



dt 

1 

1 



B(x,t) + ^D(x) 



{ax - bx 3 ) + Q(x, (Sx n )) +Acosujt 

1 



P(x,t) 



D(x) = ±[Q, 



x 2 + D a 



Eq. 



can then be rearranged into the following form 

Qix 



dP(x,t) _ ld_ 
dt 7 dx 



-ax + bx 



7 



Acosut - Q(x, (Sx n )) + ~^-(D q + Q!X 2 ) 



ld_ 

7 dx v 



P(x,t) 



(22) 



(23) 



(24) 



(25) 



The quantum nature of the stochastic dynamics is manifested in two terms; the quantum correction to the classical 
potential Q and the quantum diffusion coefficient D q characterizing the thermal bath. Ea. l|25|) is the quantum 
counterpart of Eq.(5) derived in RefSi Furthermore Ea. (|25|l formally takes care of quantum correction to all orders 
and the vacuum held effect of the bath at zero temperature. 



B. Static tilting 

In the presence of static tilting, i.e., u> = 0, the stationary solution of the Fokker-Planck equation reads as 



P (x,A) = N Q 



..2 , D i 





r fkx 2 


-)} 


exp 













The term G(x), arises due to quantum corrections to the potential and is given by the following expression 

, 7 Q(x',(Sx n )) 



G(x) = exp 



dx' 



Qj (f + % 



(26) 



(27) 



A^o is a suitable normalization constant, k = 5^7 i an< i denote the barrier height and stable minima of the 
bistable potential — § x 2 + jx 4 , respectively, given by the following expressions 



(28) 



To proceed further it is necessary to find out the quantum correction term Q(x, (5x n )) more explicitly. From Eq.Q 
we have 



Q(x, (Sx n )) 



V"'(x,t)_ {6&2) _V^ {§£3) + 



2! 



3! 



(29) 



For a bistable potential V n+1 (x,t) = for n > 4. For a minimum uncertainty state (5x 2 ) = 0(K) and therefore the 
higher order terms of the order {5x 3 ) can be safely neglected. Furthermore since the stochastic system is overdamped 
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the quantum corrections are primarily relevant in the time scale O(-) It is therefore important to calculate the leading 



order quantum correction 



V"'(x,t) 



(Sx 2 ). To this end we consider the overdamped version of operator equation (J2J) and 



use x(t) = x{t) + Sx(t), where x(t) = (x(t)) is the quantum mechanical mean value of the operator x. By construction 
[Sx, Sp] = ih and {Sx) = 0. We then make use of the overdamped version Eq. i|17fl and obtain the operator equation 
after quantum mechanical average over bath modes with product separable coherent states 



-/Sx + V"(x, t)5x + V ^-V n+1 (x, t)(Sx n - (Sx 71 )) = 



(30) 



n>2 



With the help of operator equation l|3U|) we obtain the coupled equations for (8x n ) 



j t {^ 2 ) =-- [2V"{x.t)(8x 2 ) + V"'{x,t)(Sx 3 )] 



(31) 



(Sx 3 ) = — 
at 7 



W"(x,t)(5x 3 ) + ^V"'(x,t)(Sx 4 ) - ^V"'(x,t)(Sx 2 ) 2 



(32) 



and so on. To take into account of the leading order contribution (Sx 2 ) to Q{x, (Sx n )) explicitly we retain only the 
first term in Ea. l|31|l 



d(5x 2 ) = — V"(x,t)(Sx 2 )dt 

1 



(33) 



The overdamped deterministic motion on the other hand gives jdx = —V'{x,t)dt which can be used in Ea. H33|) to 
eliminate dt and obtain after integration 



{S£2) d(5x 2 ) T V"(x,t) 



(34) 



where x c is a quantum mechanical mean position at which (Sx 2 ) becomes minimum, i. e., (Sx 2 ) Xc = ^K/ujo, wq being 
defined in Ea. l|15|) . Explicit integration yields 



where A c is given by A c 
From Ea. lfT%|) we have 



(Sx 2 ), 



\v(x c ,t)Y 



(Sx 2 ) = A c [V'(x,t)Y 



This quantum correction parameter can be estimated as follows: 



(35) 



[V'(x, t)f = x 2 [a 2 + C 2 (i) + 2a£(f) + 2(a + £(t))bx + b 2 x 
By averaging Eq. (|36|l over external noise ^ (t) using Eq. 119|l we have from Eq. 135|) 

(Sx 2 ) = A c x 2 (a 2 + Qi + b 2 x A - 2abx 2 ) 

r„2 ! o.n 1/2 

The minimum of (Sx 2 ) at be obtained from Eq. I|37|l as 



(36) 
(37) 

by neglecting b 2 term (since 

a>6). 

With the help of Ea. (|35|l it is easy to obtain from Eq.J^HJ the leading order quantum correction averaged over 
external noise as given by (we have dropped the overbar sign used in Eq. l|19|l . for convenance) 



4ab 



Q{x, (Sx n )) = -6A q bx 3 (a 2 + Qi + b 2 x 4 - 2abx 



(38) 



where A g = Q(x, (Sx n )) in Ea. l38|) is an odd function, which contributes to the potential force term, so that the 
system retains its symmetry (i.e., reflection symmetry of the bistable potential remain unchanged) 

Now we are in position to discuss the properties of the stochastic processes x(t), in terms of its steady state 
distribution function P(x,A) for different regimes of coherent and inherent system parameters. 
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In absence of the static tilting i.e., A = 0, the process x(t) diffuses on the entire x axis, the relevant distribution 
function is symmetric for x — > —x and (x(t)) = 0. As temperature of the system decreases the probability distribution 
decreases at the origin, since D q decreases with temperature. In the classical limit as T — > 0, D q — » 7&bT — * we 
observe P(x,0) — > at the origin. This is shown in solid line in Fig.l. So the process x(t) remains confined in any 
one half axis depending on the initial condition. Unlike the classical limit at T = the distribution function P(x, 0) 
in the quantum limit takes the following form; 



P o (x,0) = N 



DM 



Qi 



exp 



(1-Ai) 



(39) 



Here 7 is assumed to be unity, D q (0) is zero point diffusion coefficient having non-zero value ^huj Q , Aj is the 
contribution due to system nonlinearity correction as given by 



Ai = 6A 9 (a 2 



(40) 



Hence the process x(t) does not remain confined in the one half axis depending on the initial condition even at 
zero temperature due to quantum mechanical tunneling. This is an important aspect of bath quantization. The 

distribution function peaks at x m ~ ^ ^(l-Ax) — ~\ 2 , where n = — ^ + k (l + fc ^^ - J ■ Thus compared to 
classical limit the peaks of the distribution function shift slightly towards the origin due to zero point contribution 
of the bath oscillators. All the above observations are in the limit a > Qj and are shown in Fig.l and Fig. 2. In the 
zero temperature limit as the external noise intensity is increased to Qi = a, the distribution function Ea. l|39|) attains 



an unimodal form since 



Ay 



< 1; a phenomenon like phase transition occurs both in the classical and the quantum 



limit. This is presented in the Fig. 3. In the quantum case the distribution function deviates from the Gaussian shape 
due to quantum corrections arising out of higher order nonlinearity. It is important to note that the quantum particle 
is more localized than the classical particle over the cc-axis, a feature typical of quantum localization as shown in this 
case. For Qi > a, the distribution function Ea. H39|) tends to be singular at the origin since n < and the appearance 
of long tails is characterized by stochastic stabilization. 



C. Periodic tilting 



When the periodic tilting is present u> ^ and A ^ 0, the process x{t) is no-longer stationary and the time- 
dependent distribution function P(x, A, t) is necessary to describe its steady-state. However in the limit of low forcing 
frequency 2 ^ there is enough time for the system to reach the local equilibrium during the period of i. Then the 
quasi-steady state distribution function reads as 



Qi 



P(x) = NG(x) 
where the time r is used as fixed parameter so that 

L(x) 



exp 



kx 2 



L(x) 



A COS LOT 



(DqQl) 



arctan 



(41) 



(42) 



The quasi-steady state distribution function P(x) is plotted in Fig. 4 and Fig. 5 as a function of position coordinate 
when the strength of additive and multiplicative noises D q and Qi, respectively, are varied. The probability density 
function P{x) is depicted in Fig. 4 for several values of quantum diffusion constant and for a fixed value of A(=-0A). 



The peak located at x 



is much higher than that at x = +(f-) 2 when D q 



0.9. 



When D q is increased to 



10, the distribution function P(x) gets broadened and the left peak is decreased to almost equal height to the right 
one. Fig. 5 depicts a situation where the multiplicative noise strength is varied. It is observed that the peak of the 
distribution function at x = —(f) 5 is decreased and moves towards the origin when the strength of the external noise 
is increased from 0.08 to 0.7 with a gradual lose of Gaussian form. 



III. THE CO-OPERATIVE EFFECT; MULTIPLICATIVE STOCHASTIC RESONANCE 



In the previous section we have been analyzed the stochastic process x(t) in terms of the probability distribution 
functions. We have compared the distribution function of quantum mechanical mean values with the corresponding 
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classical one in the different parameter regime which reflect the profound effect of bath quantization at very low 
temperature. But this is not sufficient to describe completely the cooperative effect among the inherent thermal 
noise, the external noise and weak periodic signal. To analyze such nonlinear cooperative effect, i.e., the phenomenon 
of stochastic resonance in the bistable system the most common way to quantify the effect is through signal-to-noise 
ratio (SNR). SNR is defined as the ratio of the strength of 5-spikc of the power spectrum to the back ground noise 
level. Following^ the expression of SNR in the bistable system can be derived from the well-known two-state model 
approach as 




4Wo 2 



(43) 



where Wq is the Kramer's hopping rate over the activated barrier in absence of periodic force and ±xq are the two 
stable states. 

The simplest way to calculate Kramer's rate is through mean first passage time (MFPT). The exact expression for 
MPFT— for a particle to reach the final point x 2 , the quantum mechanical mean position from an initial point x\ is 
given by (W = y-) 



X 2 ) 



dx 



P s {x)dy 



I X1 D(x)P s (x)J_ c 

From Ea. (|22l - I24[l . we obtain the stationary distribution P(x) (A — and 7 = 1) as 



where 



P s (x) = ND-2(x)eyLV 



D{x) =Q lX 2 + D q ; U(x) 



U{x) f x Q(x, (5x n )) 



D(x) 



dx 



x ay - by 3 

" (W&)^ 



(44) 



(45) 



(46) 



The two points x\ and X2 are —xq and 0, respectively, 
expression for T m in the usual way 



Now using steepest-descent approximation we obtain the 



T m { Xl -» x 2 ) = 27r\U"(x 1 )U"(x 2 )\-i exp 



U( Xl ) - U(x 2 ) 
Qi 



exp [cj>{x 2 ) - <l>(xi)] 



where 



4>{x) 



x Q(y,(Sy n )) 
D{y) 



dy 



U(x) 



-~x 
2 



—x 
4 



(47) 



(48) 



Having determined the hopping rate we are in position to explain the effect of different parameters on the signal-to- 
noise ratio. From the equation the effect of increasing input signal A is quite straightforward. On increasing the input 
signal one observes the increase of SNR. The implication of varying the signal frequency is, however, complicated. 
To explore this effect we present in Fig. 6 the variation of signal-to-noise ratio as function of temperature for several 
values of signal frequency(w). The solid curve represents the value of SNR, for lowest frequency = 0.0001) and 
dashed refers to the same for the largest frequency (ui — 0.1). The SNR decreases as the frequency is increased in the 
lower temperature range. When the signal frequency is low, the particle gets time to approach a local equilibrium, 
and it has a high likelihood of hopping from upper well to the lower well during the half cycle. Hence a type of 
resonance phenomenon is expected to occur at optimum level of thermal fluctuation. At higher temperature such 
type of nonlinear cooperative effect is destroyed due to substantial probability of hopping back to the upper well. At 
the high frequency regime on the other hand such nonlinear cooperative effect is less likely to be observed, because 
a fewer number of particles find the time to hop to the lower well during each cycle. On increasing temperature one 
may increase the inter-well transition rate, but at the risk of being kicked in the anti-phase direction. Although SNR 
passes through a maximum as one varies the thermal noise strength it may be considered as some type of compromise 
between input signal and inherent thermal fluctuation rather than a resonance. 

Fig. 7 depicts the variation of SNR as a function of temperature comparing the classical and the quantum limits 
and shows that SNR is larger in quantum limit at low temperature and the difference become insignificant at high 
temperature as expected. It is also interesting to note that in contrast to classical limit SNR is not zero at absolute 
zero due to zero point contribution of the thermal bath. 
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IV. CONCLUSION 



Based on the study of stationary probability distribution function and signal-to-noise ratio we have investigated 
the problem of quantum stochastic dynamics of a system in a bistable fluctuating potential which is driven by a slow 
periodic or static field. Our conclusions are summarized as follows: 

(i) Temperature profoundly affects the characteristic of the bimodal stationary probability distribution function 
which reveals that in contrast to classical case, where the system remains confined in one of the wells, quantization 
allows a finite population at barrier top even at zero due to vacuum field effect of the thermal bath in absence of tilting 
and for a low noise strength of the modulation for fluctuating well. Thus although at zero temperature the passage 
between the wells is classically forbidden, quantum tunneling makes the passage smooth for population transfer. 

(ii) For a static tilting and a higher strength of stochastic modulation of the well, the unimodal stationary probability 
distribution becomes more localized and acquires non-Gaussian character on quantization as compared to that for the 
corresponding classical system. 

(iii) Very slow external periodic driving may localize the system in one of the wells. The peak of the distribution 
function shifts towards the origin as the strength of modulation of fluctuating potential is increased. 

(iv) The effect of quantization also makes its presence felt significantly in the signal-to-noise ratio at low temperature. 
The method presented here takes care of the statistical properties of the thermal bath in terms of Wigner canonical 

thermal distribution which remains a valid pure state distribution even at absolute zero. This help us to consider the 
vacuum field limit in the treatment of stochastic processes the governing equation of which is classical in form but 
quantum mechanical in its content. Moreover although formally one takes into account of the quantum correction 
terms due to nonlinearity of the potential to all orders in the stationary distribution function, it is possible to estimate 
these corrections order by order for actual calculation effectively to a good degree of accuracy. We thus hope that 
our observations will be useful for experimental studies on tailored quantum systems like quantum dots at very low 
temperature. 
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Figure Captions 

Fig.l: Variation of probability density function (P(x)) as a function of quantum mechanical mean position(x) in 
absence of tilting (A=0) for the parameter set a = 1.0, b = 0.06, A q = 0.02, Qi = 0.1 and (i) D q = 0.0(-fk B T = 0) and 
A q = 0.0 (solid line) [classical] (ii) D q = 0.5 (dashed line), (iii) D q = 3.0 (dotted line) and (iv) D q = 5.00 (dash-dot 
line) . 

Fig. 2: Variation of P(x) as a function of x at zero temperature comparing quantum(b) and classical(a) limits for 
the parameter set a = 1.0, b = 0.006, Qi = 0.5, A = 0.0 and A q = 0.02. 

Fig. 3: Variation of P(x) as a function of x at zero temperature limit comparing quantum(b) and classical(a) phase 
transition for the parameter set a = 1.0, b = 0.006, A = 0.0, A q = 0.02 and (i) a > Qi — 0.5 (dotted line) and 
Qi = a = 1.0 (solid line) 

Fig. 4: Variation of P(x) with x for different values of quantum quantum diffusion coefficient (D q ), D q = 0.9 (solid 
line), D q = 3.5(dotted line) and D q — 10.0 (dashed line) in presence of tilting (A=-0.1) for a = 1.0, b = 0.006, 
Q/ = 0.1, u) = 0.1, A q = 0.02 and t = 1 

Fig. 5: Variation of P(x) with x for different values of external multiplicative noise strength (Qi), Q = 0.05 (solid 
line), Qi = 0.1(dotted line), Qi = 0.3 (dashed line) and Qi = 0.7 (dash-dot line) in presence of tilting (A=-0.1) for 
a = 1.0, b = 0.006, Qi = 0.1, w = 0.1, A q = 0.02 and t = 1 

Fig. 6: Variation of Signal-to-noise(SNR) ratio with temperature(T) for the parameter set a — 2.5, b = 0.001, 
Qi = 3.0, A q = 0.02, A = 0.008 and (i) w = 0.0001 (solid line), (ii) uj = 0.001 (dotted line) and (iii) uj = 0.1 (dashed 
line) 

Fig. 7: A plot of signal-to-noise ratio vs temperature comparing quantum and classical limit for the parameter set 
a = 2.5, b = 0.001, Qi = 3.0, A q = 0.01, A = 0.008 and u) = 0.1 (i) classical case (solid line), (ii) quantum case 
(dotted line) 
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